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Abstract 

We show a text-book potential for single-field inflation, namely, the Coleman- Weinberg model 
can induce double inflation and formation of primordial black holes (PBHs), because fluctuations 
that leave the horizon near the end of first inflation are anomalously enhanced at the onset of 
second inflation when the time-dependent mode turns to a growing mode rather than a decaying 
mode. The mass of PBHs produced in this mechanism lies in several discrete ranges depending on 
the model parameters. We also calculate the effects of non-Gaussian statistics due to higher-order 
interactions on the abundance of PBHs, which turns out to be small. 
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I. INTRODUCTION 

The primordial perturbation generated in the inflationary epoch is believed to 

be the origin of large-scale structure observed in the universe today. We have accurate 
information on the primordial perturbation through the observation of anisotropy of the 



a. 



cosmic microwave background(CMB) with the help of cosmological perturbation theory 
However, the CMB observation provides us with the information on the perturbations on 
a limited range of scales. Therefore, we cannot say, a priori, anything on the primordial 
perturbation at the smaller scales from CMB data. If the perturbation at these scales is 
order unity, primordial black holes can be produced when the scale of overdensed region 
crosses the horizon |4|]. The typical mass of these black holes is given by horizon mass at 
horizon crossing^: 

Me„ ^ |p(//-')» ^ 1^. (1) 

where Mq = (SnG)-^/^ = 2.4 x 10^^ GeV is the reduced Planck mass^. Though PBH 
evaporates through the Hawking radiation process [s], those with mass greater than 10^^ g 
can remain until the present time 3, ^1 . These PBHs can be an origin of intermediate mass 

CI un 

black holes j,9| or dark matter in the Universe p, HOj. 

It is difficult to produce appreciable numbers of PBHs in the simple single-field slow-roll 
infiation models which predict nearly scale-invariant power spectrum of curvature perturba- 



tions because the amplitude of fiuctuations on small scales cannot be much differeiit 
that normalized by CMB observation on large scales. One exception among others 
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is the chaotic new infiation model [13!] where double infiation is realized with a single field. 



In Ref. 131], however, calculation of curvature fiuctuation was done using the slow- roll for- 
mula in which only the time- independent mode has been taken into account. In the present 
paper, we solve the evolution equation of each Fourier mode of fiuctuations properly and 
find an anomalous growth due to the temporal deviation from slow-roll evolution between 
two infiationary stages. This provides the first realistic example of the anomalous growth 
of perturbation in the superhorizon regime discussed in Ref. [M]. As a result we obtain a 



^ Strictly speaking, the mass of black hole depends on the density perturbations 6 and black holes with a 
small fraction of horizon mass are also produced [9| ■ Since the contributions of these black holes are small 
compared to the black holes with horizon mass, we ignore this dependence for simplicity. 

^ We use units for which c = 1. 
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power spectrum highly peaked on some scales depending on the values of the model param- 



eters and PBH formation is more easily realized than concluded in Ref. [13|. We search for 
the values of parameters with which appreciable numbers of PBHs are produced under the 
observational constrains [isl. We also analyze the effects of non-Gaussianity generated in 
this model on the abundance of PBHs. 

The organization of this paper is as follows. In section [Tll we introduce the chaotic 
new inflation model and explain the background evolution in this model. In section llllt a 
mechanism of the enhancement of the perturbation is explained and the power spectrum 
in the chaotic new inflation model is given. In section IIVI we give an expression of PBH 
abundance resulted from a peaked power spectrum. In section |Vl we estimate non-Gaussian 
correction to PBH abundance. In section Wl\ we calculate the PBH abundance with various 
values of parameters and give a relation between mass and PBH abundance. Section IVIII 
contains conclusion of this paper. In this paper, we use curvature perturbation in the 
comoving gauge ( as a degree of freedom of scalar perturbation. 



II. CHAOTIC NEW INFLATION MODEL 



We consider a single-field inflation model with the Coleman- Weinberg potential 
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Historically, this potential was first used to realize new inflation 17|], but phase space 
consideration has led to the conclusion that chaotic inflation ISj is much more likely to 



occur 



19| . So we start with a large field value. 



In this model, inflation can occur twice |20j. First, chaotic inflation occurs. After chaotic 
inflation, inflaton oscillates between the two minima of the potential. If the parameter v is 
appropriately chosen, the inflaton moves slowly in the neighborhood of the origin after the 
oscillation and new inflation occurs. For example, the number of e-folds of new inflationay 
expansion, A^new, with \H\ < H"^ is larger than 10 for v = 1.103Mg - 1.132Mg, and it 
satisfies A^ncw ^ 60 for v = 1.114Mg — 1.122Mg- According to the number of oscillation 
cycles, the parameters where new inflation occurs are distributed at certain intervals. In the 
above example, if settles to the positive potential minimum Lp = v ii v > LllOMc and to 
the negative potential minimum ip = —v if f < 1.118Mg without oscillation. New inflation 
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with A^new > 10 occurs after a half cycles of oscillation for v = 0.35510Mg — 0.35524Mg! 
and after a cycle of oscillation for v = 0.266665Mg' — 0.26670Mg' and so on. We show 
the evolution of Hubble parameter H and the inflaton (p with the values of the parameters 
(A, v) = (5.4 X 10"^^, 0.355139Mg) in Fig. [H which shows that the inflaton moves slowly 
in the neighborhood of the origin, and new inflation occurs. With these values of the 
parameters, (p settles to the negative potential minimum after an oscillation. 
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FIG. 1: The evolution of Hubble parameter(left) and the inflaton(right) with the values of the 
parameters (A, v) = (5.4 x 10"^^, O.SSSlSQMc). is a value of scale factor at the initial time. 

In describing the evolution of the inflaton during inflation, it is convenient to introduce 
the following Hubble slow-roll parameters: 



H 



V 



e 

We' 



(3) 
(4) 



where dots denote differentiation with respect to the cosmic time t. We show the evolution 
of these slow-roll parameters in Fig. O The values of the model parameters are the same as 
those employed in Fig. [H We can see that slow-roll conditions are not satisfied while the 
inflaton is oscillating between the two minima. As we see in the next section, the existence 
of this period is important for enhancement of curvature perturbation. 

III. ENHANCEMENT OF CURVATURE PERTURBATION 



Curvature perturbation is enhanced in the models where slow-roll conditions are tem- 
porarily broken as stated in Ref. jl^. In this section, we briefly describe the mechanism of 
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FIG. 2: The evolution of slow-roll parameters with the values of the parameters (A, v) = (5.4 x 
10-l^ 0.355139Mg). 

the enhancement and give a power spectrum of curvature perturbation in the chaotic new 
inflation model. 



A. Evolution of curvature perturbation 

Curvature perturbation in the comoving gauge (, in terms of which the amplitude of 
perturbation in the intrinsic spatial curvature of the comoving slicing TZc is written as 

= ^V\, (5) 

evolves according to an equation 



''^' + (3- + .)^H-f^V^.-0, (6) 



where N is the number of e-folds and Cfc is the Fourier transform of (: 



a = jd'xat,x)e~''''. (7) 

In the slow-roll inflation regime, the coefficient of the second term in Eq. ([6]) is positive 
and therefore the solutions of Eq. in the long-wavelength regime, where the last term 
is negligible, are a constant mode and a decaying mode. The time derivative d(k/dN, 
which corresponds to the decaying mode, diminishes in proportional to a~^. If we neglect 
the last term in Eq. ([6]), d(k/dN diminishes as a~^. However, the last term decreases as 
[k/aHYQ]^ oc a^^. As a result the second term in Eq. ([6]) soon becomes comparable to the 



last term in Eq. ([6]), and d(k/dN diminishes as a~^. Therefore Cfc soon becomes constant 
after horizon crossing and the power spectrum of curvature perturbation is given by the 
squared amphtude of the vacuum fluctuation^ 

1 k-^ 



at the time the mode crossed the horizon: 



2eM2 a2 



(8) 



k'Kk 



27r2 
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(9) 



k=aH 



The power spectrum given by Eq. ([9]) is well approximated by the form V(^{k) oc fc" ^ 
where n — 1 = —2e — rj, namely, a nearly scale-invariant spectrum with small value of 



slow-roll parameters. Since the observec 
generically gives spectral index n < 1 



amplitude is small and single-field inflation models 



2l| . this power spectrum leads to a small number 



of primordial black holes. On the other hand, the coefficient of the second term in Eq. 
can be negative in the model where slow-roll conditions are temporarily broken such as the 
chaotic new inflation model, in which we find e ~ 3 and 7] < near the end of the oscillatory 
phase so that 3 — e + ?7<0. In this case, Eq. (jS]) has a growing mode solution instead 
of the decaying mode solution outside the horizon. Therefore Cfc grows even after horizon 
crossing and its amplitude is enhanced. This enhancement of the perturbation enables even 
a single-field infiation model to produce a large number of PBHs. 



B. Power spectrum of curvature perturbation in the chaotic new inflation model 

We have estimated the power spectrum of curvature perturbation in the chaotic new 
infiation model by solving Eq. numerically. In Fig. [3l we show the power spectrum 
which is normalized to the amplitude observed by WMAP [3I, 2 x 10^^ at A; = 0.002/Mpc 
by choosing A appropriately. The values of the model parameters are the same as those 
employed in Fig. [T] and Fig. [2J 

^ In subhorizon scales, is given by the quantum fluctuation in the vacuum state. Therefore, Cfe^'Cfc 
the minimum value 1 under the uncertainty relation CkT^Ch — 1- Here, tt.^^, is the momentum conjugate to 
Cfe and expressed as 2a^eMQ)^^. Then, approximating ^ {k/a)(k, we obtain the amplitude ([5]). 
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FIG. 3: Power spectrum of curvature perturbation (solid line). This spectrum is calculated under 
the parameters (A, v) = (5.4 x 10~^^, 0.355139Mg). We show also a power spectrum estimated 
by using the formula ([9]), which is used for a slow-roll inflation model (dashed line). 

In Fig. [31 it is observed that the power spectrum deviates from the one estimated by 
using the slow-roll formula iQ and the enhancement of the perturbation has occurred. The 
power spectrum has a peak with amplitude ~ 6.2 x 10^^ and results in formation of a large 
number of PBHs with mass corresponding to the scale of the peak. The scale of the peak 
corresponds to the scale which crossed the horizon near the end of chaotic inflation. The 
non-constant mode for the scales which crossed the horizon earlier is a decaying mode during 
slow-roll inflation, and becomes exponentially small as a^^. Therefore, even if it turns to 
a growing mode temporarily after the flrst inflation, only the modes which left the horizon 
in the late stage of the flrst inflation are enhanced to a visible level. The scale and the 
amplitude of the peak vary with the parameter v. For increasing number of the oscillation 
cycles of inflaton, a larger amplitude is obtained. We investigate the abundance of PBHs 
resulted from the peaked spectrum in the following sections. 



7 



IV. ABUNDANCE OF PRIMORDIAL BLACK HOLES 



In this section, we give an expression of the PBH abundance resulting from a strongly 
peaked power spectrum. 

We estimate the PBH abundance based on the Press-Schechter method 



22 



23|. In 



this method, a PBH with mass greater than M is formed when the perturbation which 
is smoothed on scale Rm corresponding to M exceeds the threshold Cth- The smoothed 
perturbation C,r^j is defined by 

Cr^^{x) = j dV W{\x' - x\/RM)ax'). (10) 

where W{x/R) is a window function. The fraction of the energy density of the Universe 
collapsing into PBHs with mass M < Mbh < M + AM at the time they form is given by 

P(M: AM) . = _2 dM / dC (11) 



Ptot Jm icth 

where the prefactor 2 is due to Press-Schechter's prescription and Pr^,^ is the probability 
distribution of ^Rm independent of x because of homogeneity of the universe. 

Therefore, we omit the argument x in the following. For the moment, we assume to be 
Gaussian distributed, which is the case to the lowest order of perturbation, and consider 
non-Gaussian correction in the next section to test the validity of this assumption. Then 
Prj^ is Gaussian with variance 

- / fwikRrV,, (12) 

where W{kR) is the volume-normalized Fourier transform of the window function W{x/R). 

We estimate the PBH abundance, ffTTl) . resulting from a strongly peaked power spectrum 
at the scale = fcpg^k- Since the power spectrum is strongly peaked we can approximate 
it to be 5-function in Eq. ( |T2|) . This approximation is adequate, for Pr^^ is exponentially 
sensitive to cx/j^^ and a small difference in ct/j^j- leads to a large difference in Pr^^ so that only 
the value of aj^ at the peak is important. Under this approximation and using a top-hat 
function as W{kR), we obtain daj^^JdM ~ V^i^kpcakj^iM — Mpeak) where Mpeak is the mass 
corresponding to the scale of the peak. Therefore, only PBHs whose mass is Mpeak form. In 
this approximation, its abundance is given by 

/3(Mpeak) = ^ ^ dC (C' - l)e-^'/^ (^C = C/v/^c(W)) , (13) 



where /3 is independent of AM under the current approximation. In the peaked spectrum, 
the amphtude at the scale of the peak can be large while the amplitude at the observed scale 
is consistent with the observational value. Therefore, a large number of PBHs are produced 
from strongly peaked power spectrum. 



V. NON-GAUSSIAN CORRECTION TO PBH ABUNDANCE 



The amplitude of perturbation producing PBHs is so large that non-Gaussinity of cur- 
vature perturbation due to higher-order effects can be important. In the case of slow- roll 
inflation, this effect on the formation of PBHs has been studied in Ref. [24] and shown 
to be negligibly small contrary to the intuitive expectation. However, in models where 
slow-roll conditions are temporarily violated such as the chaotic new inflation model, large 
non-Gaussinity arises 25j and can modify the PBH abundance obtained above. In this 
section, we investigate this possibility. 



A. Three-point correlation functions 

Here in order to estimate the effects of the deviation from Gaussian due to the higher- 
order interaction in the action, we first evaluate three-point correlation functions. 

We sketch the derivation of the three-point correlation functions j26|. The three-point 
functions are estimated by perturbative expansion with respect to the interaction. For the 
estimation of the three-point correlation functions, the cubic terms of C ^ in the action give 
relevant interactions, which consist of 

^3 = Jdt Ls{CX;t) 

= Ml Utd?x [aeC(VC)^ + a^^ei^-^C^ - 3a\CC^ 

J (14) 



where if) is defined by 



^Ij^-^ + aW-'C- (15) 



^ To be precise, this C is a generalization of C used in the linear perturbation theory [26] . At linear order, 
two Q coincide with each other. 
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From these terms, the interaction Hamiltonian up to the third order reads 

= -L3{CiXi;t) 



— (30 - iJ-^C/)(V.V,-V'7V,V,^, - V^V'/VV/) - 2a-V,V'/V,C/V2^,] 



(2vr) 

+ n Hdt Hdt^"''^ Hdt Hdt Hdt 



(16) 



where vr^ is the momentum conjugate to ( and the variables with subscript / denote variables 
in the interaction picture. Coefficients Ti*^*-* (i = 1,2,3,4) are given by 







^(2) 



-eA;^(sin^ 9/z + cos ^) + y sin^ 6, 



7^(3) _ 3g ^ ^(sin^ ^ + 2^ cos 9) + ^k\z - 2 cos0) sin^ 9, 

7^(^) = -e--sin2^, 
2 



(17) 
(18) 
(19) 
(20) 



where we have defined 



k = cos9 = {ki ■ k2)/kik2, z = kl/kik2, 

ari 



(21) 



for simplicity of expression. For an equilateral triangle, ki = k2 = k^ = k, the values of 
these quantities are 

k = —, cos6' = --, z = l. 
aH 1 



In calculating the three-point correlation functions, the in-in formalism 27|] is used since 
we want to calculate an expectation value with respect to the 



(22) 



where Sq is a positive infinitesimal constant^. To leading order of Hj, the three-point corre- 



In an abuse of language, we use the same symbol Cfc for denoting a quantized variable as a classical one. 
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lation function is given by 



to 

= 2 / dt'Im i{aAt)CkAt)Cksiit)HjiCki, vrc,/; t'))) , (23) 

J to 

where the t' integration contour is deformed so that both bra and ket are projected on to the 
vacuum state at t ^ — oo. According to the behavior of (kj, the integrals in Eq. fl25]) can be 
sphtted into three parts, an integral over the region inside the horizon, the region around 
horizon crossing and the region outside the horizon. In the former two parts the deviation 
from the slow-roll inflation models is small, since the third term in Eq. ([6]) is dominant in 
these parts. The deviation arises in the last part, because Cfc grows outside the horizon in 
the chaotic new inflation model and not in the slow-roll inflation models. In the slow-roll 
inflation models, the contribution from the part outside the horizon is negligible since each 
term in the interaction (fT6!) includes (i or k/aH which have small values outside the horizon 
and, furthermore, the commutator of (^s or its time derivatives vanish as a"*^ (z/ > 2) ^. 
The three-point correlation functions are suppressed by the slow-roll parameters estimated at 
horizon crossing In contrast, in the chaotic new inflation model, (j has a non-negligible 
value outside the horizon and the integrals over the region outside the horizon contributes 
to the three-point correlation functions. The reason for obtaining large non-Gaussianity is 
different from that in Ref. [3]. In the chaotic new inflation model it is the growth of the 
perturbation outside the horizon, while in Ref. 25(] it is the characteristic behavior of the 
slow roll parameters near or inside the horizon. 



B. Correction to PBH abundance from three-point correlation functions 

In the following, we give an expression of the PBH abundance resulting from the pertur- 
bation with the three-point correlation functions. 

The probabihty distribution Pr^^ can be expressed as 

PrACrm) = ^J dVR,, $(r7R,Je-''-MC«M, (24) 

^ If Cki{t) is constant Cki{t) = CkF^^j the three-point functions vanishes because a term 

Ck,imk,iit)Ck,i{t)a,iit')CkAt')Cki(t') becomes a real number laTCkT'C^Tl' ^^is case. Then 
the leading term of the commutator of Cfc's have a decaying mode in Ck's, and decreases as a~^. With 
time derivatives of Ck's, the commutators decreases more rapidly. 
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where $ is defined by 

^r]Rj = (e^^«MC«M). (25) 
$ can be expanded by cumulants of (rm^ {(rm^^^- 

'^(VRm) = exp ^^^(a)c) . (26) 

\m=0 ■ / 

Tlie cumulants of Cr^^ can be expressed by tlie connected part of the correlation functions 
of (k- With the formula I^E^, we can estimate the correction to Pr^^ from higher-order 
correlation functions. 

The correction from the three-point correlation functions can be calculated, retaining 
terms up to m = 3 in Eq. (!26l) . In the case |(C^^f)c| is much smaller than we can give 
a concrete expression of the corrected PBH abundance. In this case, we can expand $ with 
respect to J which is defined by 

j^USkill^ (27) 



and get 



Km 



6 ctr 



1+ ( 3C^_ ^ I J 



e^P ( ) • (28) 

Rm ■ 



J can be expressed by the three-point correlation functions of Cfc as 

J = ^ I (^'h I d^A;2 I d^fca W{hR)W{hR)W{hR){CuXM ■ (29) 

R]\£ J J J 

Substituting the probability distribution Pr^^ given by Eq. fl28l) to Eq. f|TT]) . we get the PBH 
abundance including the correction from the three-point correlation functions. 

Because of homogeneity of the universe, the three-point correlation functions of Cfc can 
be written as 

{(ikAt)CkMCk,{t))c = Aiki, k2, k3)5\ki + k2 + ks). (30) 

The isotropy of the universe guarantees that A{ki, k2, k^) is in fact a function of "shape" 
and "size" of triangle spanned by fej = 1, 2, 3) such as the quantities (l2TI) and independent 
of direction of the triangle. In the chaotic new inflation model, Cfc grows most significantly 
at k = /cpeak and A{ki, k2, k^) has the largest value at ki = k2 = = fcpeak- Therefore, as 
before, we can take approximation d(C|^^)c/dM ~ 47rfc|^g^]^Apeak5(M - Mpeak), where Apeak 
represents A{ki, k2, k^) estimated at ki = k2 = ks = kpeak- Under this approximation, the 
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PBH abundance is given by 

1 



peak J 



Cth 



(C'-i)-(C'-8C^ + 9C)J, 



peak 



C = C/Vn(Vak) 



J neak — 



2vr ^poak^peak 



where 

The PBH abundance with some values of Jpeak is shown in Fig. HI 

log a 




log Pg 



(31) 



(32) 



FIG. 4: PBH abundance with some values of Jpeak- The values of Jpeak are —0.1, —0.01, 0, 0.01, 0.1 
from below. 



C. Estimation of correction to PBH abundance in chaotic new inflation model 

With the formulae obtained in the previous subsections, we can estimate the correction to 
the PBH abundance from the three-point correlation functions in the chaotic new inflation 
model. We have carried out the integration in Eq. (1231) numerically over the region outside 
the horizon, where the deviation from the slow-roll inflation models arises. In the other 
regions the contributions are almost the same as those in the slow-roll inflation models and 
do not modify the PBH abundance relevantly. Since the calculation is done only outside the 
horizon, we do not have the problem in implementing the ie prescription numerically. 
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For the values of the parameters (A, v) = (5.4 x 10~^^, 0.355139Mg) with which the 
amphtude at the peak has the values for producing relevant number of PBHs, ~ 6.2 x 10~^, 
the value of Apeak is estimated to be 

Ck^peak ~ 10-^ (33) 

On the other hand, the value obtained by estimating the contribution from the region around 



horizon crossing as done in slow-roll inflation models is calculated to be 



fc^e^k^peak = 487r7NLPc ~ 10" ' (34) 

where the value of V(; is estimated at horizon crossing, ~ 10^^^ (see Fig. [3]). /nl is an 
estimator usually used for parameterizing the size of non-Gaussianity observed in the CMB 



261] and has the value of the order of the slow-roll parameters estimated at horizon crossing, 
~ 10^1 (see Fig.E]). 

We get larger three-point correlation functions than those without the enhancement of 
the perturbation. However, the parameter Jpcak which is an estimator for the correction to 
the PBH abundance from the three-point correlation functions is estimated to be 

Jpeak ~ 10-^ (35) 

This value is too small to modify the PBH abundance relevantly. Though large three-point 
correlation functions are obtained, since a denominator of the estimator is also large, only 
a small correction is obtained. 

The correction to the PBH abundance from A^-point correlation functions appears with a 
factor V^^^'^ . Therefore the corrections are expected to be small as the correction from 
three point functions are. We can therefore use the expression (fT3l) for the PBH abundance 
safely. 



VI. PARAMETER SEARCH 

In §1111 we have shown the curvature perturbation is enhanced in the chaotic new 
inflation model. In this section, we calculate the PBH abundance using Eq. (fT3ll with 
various values of the parameter v and search the parameter with which relevant number of 
PBHs can be produced. For each f , the parameter A is fixed by the power spectrum of C 
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observed by WMAP, = 2 x 10 ^ at k = 0.002/Mpc. In the linear perturbation theory, 
the power spectrum V(^{k) scales as^ 

Vcik) ^ V'^{k) = XVc(X~^^^k) for V ^ V = XV. (36) 

Since the horizon mass at the matter-radiation equality is ~ 10^'^ Mq, PBHs that can be 
the origin of intermediate mass black holes or dark matter are produced in the radiation- 
dominated epoch. In the radiation-dominated epoch, the threshold value of the density 
perturbation 5th is given by 1/3 0,Q]. Then, we find the corresponding value of the curvature 
perturbation (^th = 0.75 with the help of the formula in the linear perturbation theory, 
5k=nH = 4(^fc=a///9 [23]. A similar value has been obtained by numerical calculation in 



Ref. m 



. Equation ([T]) relates the mass of PBH produced in the radiation-dominated epoch 



to the scale of the perturbation as 



10^GeV\ Va*^A/ 



^ = 5.4g (^TTvF) (^) , (37) 



where the variables with subscript "*" denote those estimated at arbitrary time in the 
radiation-dominated epoch. We assume the period of reheating is negligible and have chosen 
the time at which inflation finished as time at which the variables with subscript "*" are 
estimated. We write these variables with subscript "end" . 

We have found that a relevant number of PBHs can be produced around v ~ 0.266Q7Mg, 
V - 0.35514Mg and v - 0.35522Mg. In these values of v, 1/^^ ~ {Xv^y/^ ~ 10^^ GeV and 



Eq. (1371) gives 

M~10^g fl^l , (38) 

V -Kend / 

where -Rend = l/(flend-f^cnd), which is the comoving scale crossing the horizon at the end of the 
entire inflation. The PBH abundance (3 is plotted as a function of the mass corresponding 
to the scale of the peak Mpeak in Figs. [5] and [6] with the observational constraints. The 
constraint on PBH abundance with mass above 10^^ g is obtained from the condition that 
the abundance of PBHs is less than that of matter today. The other constraints are obtained 
from the consistency with nucleosynthesis and 7-ray observation. In Fig. [6l we can observe 



^ The equations in the hnear perturbation theory and the initial conditions for V' — XV can be written 
as those for V by defining {N' , H' , tp' , k' ,C) = iN,X-^/^H,ip,X-^/^k,\^/'^C)- Then, the relation §^ is 
obtained. 
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that PBHs with mass 6 x 10^° g can constitute a large part of matter in the universe. Though 
PBHs with some range of mass are observationally excluded to be the dominant component 
of dark matter 29|], there are no strong constraints in mass range 10^° g — 10^^ g {10~^^Mq — 
IO^^Mq). The mass 6 x 10^° g are in this range. Therefore, we can obtain PBHs that can 
be an origin of dark matter in the chaotic new inflation model. 
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FIG. 5: PBH abundance /3 as a function of mass corresponding to the scale of the peak Mpeak asso- 
ciated with the values of parameter v = 0.2666694Mg — 0.2666745^^^ (solid line). The constraints 
on PBH abundance are also depicted. 
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FIG. 6: PBH abundance /? as a function of mass corresponding to the scale of the peak Mpeak 
associated with the values of parameter v = 0.355136Mg — 0.355141Mg (solid line) and v = 
0.355217MGr — 0.355223Mg (dashed line). The constraint on PBH abundance is also depicted. 



Finally, we estimate scalar spectral index n and tensor-to-scalar ratio r and check the 
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consistency with the observational data |3|: 

n = 0.986 ± 0.022, r < 0.43 (95%CL) at k = 0.002/Mpc. (39) 

In each parameter regions, scalar spectral index n and tensor-to-scalar ratio r is estimated 
to be 

n = 0.939 - 0.931, r = 0.33 - 0.37 for = 0.266669Mg - 0.266674Mg, 
n = 0.928 - 0.922, r = 0.38 - 0.42 for = 0.355136Mg - 0.355141Mg, 
n = 0.931 - 0.924, r = 0.37 - 0.41 for = 0.355223Mg - 0.355217Mg, 

especially 

n = 0.925, r = 0.40 for v = 0.355140Mg (Mpeak = 6 x 10^° g). 

At the first parameter region, n and r are consistent with 3-year WMAP data at the confi- 
dence level of 99.9%, but lie out of the preferred range of 5-year data. At the other regions, 
n and r also lie out of the preferred range of 5-year data. 

VII. DISCUSSION 

In this paper, we have shown that the enhancement of the perturbation occurs and a large 
number of PBHs are produced in the chaotic new inflation model. A growing mode plays 
an important role in the enhancement. Due to the growing mode, there is a contribution 
to the three-point functions from the region outside the horizon. We have estimated this 
contribution and its effect on the PBH abundance. As a result. We have found that non- 
Gaussian correction to the PBH abundance due to higher-order effect is small. Further, 
we have calculated the PBH abundance with various values of the parameter v. We have 
obtained the relevant number of PBHs around v ~ 0.26667Mg, v ~ 0.35514Mg and v ~ 
0.35522Mg. In the second parameter region, the produced PBHs can constitute a large part 
of dark matter in the Universe. In the first parameter region, though the produced PBHs 
cannot be observed, we can observe the enhancement of the perturbation with gravitational 
waves generated from scalar perturbations through non-linear couplings [3^. The scale of 
the peak which gives mass M = 1.4 x 10^'^ g corresponds to GW frequency^ ~ 100 Hz, 



The power spectrum of the induced GW have the peak at kpeak/V^n [30 1 
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and power spectrum for the induced GW is estimated to be^ ~ V'^ ~ 10 ^. Therefore, the 



32|, LIGO y, TAMA |3|, 



induced GW can be detected by detectors such as GEO600 
VIRGO 0- 

In both parameter regions, however, the values of spectral index and tensor-to-scalar ratio 
are not in the preferred region of the observational data. This is because, in the chaotic 
new inflation model, the observed perturbations are produced at the chaotic inflationary 
epoch with a potential close to the quartic one, which is not preferred by the observational 
data. In this paper we adopted the Coleman- Weinberg potential as a simple model to drive 
inflationary dynamics without specifying particle physics background. In reality, we expect 
supergravity corrections to the scalar potential is important, especially for Lp > Mq- If such 
corrections modifies the potential in an appropriate way, this scenario could be a feasible 
one. Another way to rescue it is to adopt a smaller self coupling A using the scaling law fl36l) . 
so that it predicts smaller amplitudes of both scalar and tensor fluctuations on large scales 
observed by CMB etc., assuming that the observed scalar perturbations on these scales were 
created through a different mechanism, say, curvaton [36], or modulated rehating js?]. Since 
the peak we have found is so prominent that one could produce observable amount of PBHs 
even in such a case with smaller A. 
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